First-Matsubara-frequency rule in a Fermi liquid. Part I: Fermionic self-energy 
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We analyze in detail the fermionic self-energy E(a;, T) in a Fermi liquid (FL) at finite temperature 
T and frequency uj. We consider both canonical FLs - systems in spatial dimension D > 2, where 
the leading term in the fermionic self-energy is analytic [the retarded lmS^{uj,T) = C{u^ +7r^r^)], 
and non-canonical FLs in 1 < _D < 2, where the leading term in ImE^(ct;, T) scales as or ui^ . We 
relate the uP + tt^T^ form to a special property of the self-energy - "the first-Matsubara-frequency 
rule", which stipulates that Y,^{i-KT,T) in a canonical FL contains an 0{T) but no term. We 
show that in any D > 1 the next term after 0{T) in T.^{i-KT, T) is of order T'^ (T^ InT in D = 3). 
This T° term comes from only forward- and backward scattering, and is expressed in terms of 
fully renormalized amplitudes for these processes. The overall prefactor of the term vanishes 
in the "local approximation" , when the interaction can be approximated by its value for the initial 
and final fermionic states right on the Fermi surface. The local approximation is justified near a 
Pomeranchuk instability, even if the vertex corrections are non-negligible. We show that the strength 
of the first-Matsubara-frequency rule is amplified in the local approximation, where it states that 
not only the term vanishes but also that T,^{i-KT, T) does not contain any terms beyond 0{T). 
This rule imposes two constraints on the scaling form of the self-energy: upon replacing oi by iirT, 
ImE^(cj,r) must vanish and ReS^(a;,r) must reduce to 0{T). These two constraints should be 
taken into consideration in extracting scaling forms of E^(aj,T) from experimental and numerical 
data. 

PACS numbers: 71.10 Ay, 71.10. Pm 



I. INTRODUCTION 

Properties of single particle and collective excitations 
in strongly interacting electron systems continue to at- 
tract substantial interest of the condensed-matter com- 
munity. This interest is stimulated by the avalanche of 
discoveries of new materials, many of which fall into a 
category of strongly correlated electron systems, and by 
advances in experimental techniques, which allow one 
to extract, with good accuracy, the single-particle self- 
energy from angle-resolved photoemission (ARPES) data 
and the two-particle or "optical" self-energy from the real 
and imaginary parts of the optical conductivity. 

One of the most actively explored directions in the 
study of strongly correlated electron systems is a search 
for non-Fermi liquids (non-FLs) - systems in which elec- 
trons interact so strongly that they completely lose co- 
herence. Many newly discovered systems were classified 
as non-FLs because their electron spectral functions, ex- 
tracted from ARPES, are quite broad. However, a broad 
spectral function is an indication, but not the proof, that 
the system in question is a non-FL, as the Landau crite- 
rion for the FL only requires that the spectral function 
must be sharp for fermions in the immediate vicinity of 
the Fermi surface (FS). A mathematical formulation of 
this requirement is that the imaginary part of the re- 
tarded self-energy ImS''*(ti;) must be much smaller than 
w -I- 'ReY.^{uj) at the smallest w. This does not preclude 
that at higher frequencies ImE^(Lij) can become compa- 
rable to + Rel]^(a;) or even exceed it. 

To satisfy the Landau criterion, ImE^(a;) has to scale 
as oj^"*"" with a > 0. The original argument by Landau, 



based on the Pauli principle and the assumption of an- 
alyticity, yields Iml]^(cj) oc uP , i.e., a = 1. Microscopic 
calculations show that Iml]^(w) does indeed scale as uP 
in a 3D FL. The same holds for all "fractional" dimen- 
sions D > 2. For D < 2, the analyticity is, however, 
broken: ImE^(u;) scales as a;^ln|a;| in D = 2 and as 
\uj\^ in D < 2. Still, by Landau criterion, these systems 
are FLs, as long as D > 1. Hereafter we refer to sys- 
tems in which ImY,^ (uj) cx oj^ as "canonical FLs" , and 
to systems in which ImE^(cj) cx with < a < 1 as 
"non-canonical FLs" . 

The goal of this paper is to analyze the form of the self- 
energy in both conventional and non-conventional FLs at 
finite frequency uj and temperature T. We will be par- 
ticularly interested in how general is a certain property 
of the self-energy, which we will be referring to as the 
"first-Matsubara-frequency rule" or, for brevity, as the 
"first-Matsubara rule". This rule states that the self- 
energy E(wm, T), evaluated at discrete Matsubara points 
'-^m = 7rT(2TO -|- 1), exhibits a special behavior at the 
first fermionic Matsubara frequency, ujq = ttT, namely, 
E(7rr, T) does not contain terms higher than T. (The 
same happens at tj_i = — vrT.) This rule was proven in 
the past for particular cases of the electron-phonon^'^ and 
screened Coulomb^'"' interactions. In the former case, 
this rule is sometimes being referred to as a "Fowler- 
Prange theorem".^ 

Although the first-Matsubara rule operates on the 
imaginary frequency axis, it is relevant to properties of 
physical fermions with real frequencies: it requires that 
the retarded self-energy E^(cj,T), with w replaced by 
ittT, should not contain terms beyond 0{T), and thus 
imposes a constraint on the interplay between the lo and 
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T terms in 5]^(w,T). 

A 3D FL provides a simple example of how the first- 
Matsubara rule works. To order oJ^ and T^, we have 
in this case ReE^(ci;,T) = Aw, with no wT term, and 
Iml]^(w, T) = Ciu"^ + TT^T^), with a factor of exactly tt^ 
in front of T^. Ai lu — inT, ImE^ vanishes and ReE^ 
becomes of order T, hence the total T,^{iTTT,T) contains 
only an 0{T) term but no term. 

In this paper, we analyze the validity of the first- 
Matsubara rule beyond the conventional FL paradigm. 
The proof of this rule in prior work^^* was based on 
demonstrating the nullification of the leading term in 
imaginary part of the self-energy at the first Matsub- 
ara frequency. We show here that the first-Matsubara 
rule does not hold beyond the leading order for con- 
ventional FLs, and does not hold at all for unconven- 
tional FLs. Our primary finding is that E^(i7rT, T) 
scales as in all D (with an extra InT factor in 
= 3); however, the consequences of this finding are 
different for conventional and unconventional FLs. For 
conventional FLs, i.e., for 2 < D < 3, the term 
is still subleading to T^, and thus the first-Matsubara 
rule holds to order T^. For unconventional FLs, i.e., 
for 1 < D < 2, the term is of the same order 
as the leading terms in ImE^(cLi,T), and thus the the 
first-Matsubara rule is violated. In _D — 2, which is a 
marginal case between conventional and unconventional 
FLs, ImE^(a;,r) cx {uj^ + Tr^T^)\n\uj\+0{oj'^ ,T'^). While 
the logarithmic term vanishes at = iirT, the term 
does not. As a result, the first-Matsubara rule is satisfied 
to logarithmic accuracy but not beyond. 

We find that for 1 < £> < 3, the T° term in T.^^iinT, T) 
is universal, i.e., independent of the upper cutoff of the 
theory. Furthermore, its prefactor is expressed via ex- 
act spin and charge components of the forward- and 
backscattering amplitudes. 

At the same time, we find that the first-Matsubara 
rule holds to all orders in T in both conventional and 
non-conventional FLs, if the effective interaction between 
fermions, which includes dynamic screening by particle- 
hole bubbles, is assumed to connect only the states right 
on the Fermi surface. Hereafter we refer to this approxi- 
mation as the "local approximation" , as it is generally 
valid when bosons which mediate interaction between 
fermions are slow compared to fermions 

We show that, within the local approximation, the 
first-Matsubara rule relies only on the analytic proper- 
ties of the local susceptibility. For the electron-phonon 
interaction, this approximation is a key ingredient of the 
Eliashberg theoryjiS and the small parameter which con- 
trols this approximation is the ratio of the Debye fre- 
quency to Fermi energy. We consider here the case of 
an electron-electron interaction. In certain limits it can 
be approximated by an effective interaction mediated by 
collective modes of fermions in the spin or charge chan- 
nel. The collective modes are generally not slow com- 
pared to fermions themselves (their velocity is of order 
of the Fermi velocity), but they do become slow near a 



Pomeranchuk instability, when the correlation length for 
critical collective modes diverges. As a result of this di- 
vergence, the system generates a low-energy scale, below 
which near-critical collective modes become overdamped 
and slow down.^'l- The local approximation for collective 
modes is a necessary but not sufficient condition for the 
Eliashberg theory, as the latter also requires vertex cor- 
rections to be small. In the case of collective modes, 
vertex corrections are not controlled by the same pa- 
rameter which makes the local approximation valid)^^— 
and are not necessary smallJ^ii^ii^"— We show that the 
smallness of vertex corrections is not required for the 
first-Matsubara rule to work - the local approximation is 
sufficient. 

We analyze the local approximation in more detail and 
show that the first-Matsubara rule imposes two condi- 
tions: I) ImE^(z7rT, T) vanishes to all orders in T, and 
2) ReE^(i7rT,r) contains an 0{T) term but all higher 
order terms in T vanish. These two conditions are non- 
trivial because, beyond the conventional FL paradigm, 
E^(a;,r) cannot be obtained from the T = result by 
a simple replacement uj — ?► Vcj^ + tPT^. This is true for 
the subleading uj'^ terms in a 3D FL, and also for the 
leading w^, terms in non-conventional FLs. In par- 
ticular, ImE^(a;,T) in non-conventional FLs has a com- 
plex form which is very different from (w^ + tt^T^)^/^, 
and ReE''*(a;,T) also contains a complex dependence on 
ijj and T at order w^, in addition to the Aw term. Never- 
theless, as long as the local approximation is applicable, 
ImE^(w, T) vanishes at w = iirT, and ReE^(i7rT, T) re- 
duces to inXT. 

Finally, we show that the first-Matsubara rule holds 
within the local approximation even for a non-FL, e.g., 
for a system in 13 < 3 right at a Pomeranchuk in- 
stability, except that in this case the coefficient A in 
E^(i7rr,T) = inTX diverges as T -> 0. In particular, 
the first-Matsubara rule holds for a marginal FL and for 
an itinerant 2D system at a nematic quantum critical 
point (QCP). 

The rest of the paper is organized as follows. In Sec.llll 
we review the derivation of the single-particle self-energy 
to order and w^ in a conventional FL, and show where 
the relation between the and terms comes from. 
In Sec. Illli we discuss the self-energy outside of the con- 
ventional FL paradigm. We show that, in general, the 
self-energy contains terms of order T^, which do not 
vanish when w is replaced by iirT. The case of 13 = 2 is 
marginal, and we consider it separately. In Sec. IIV[ we 
discuss the self-energy within the local approximation. 
We show that, at order T^, 'S^{uj,T) has quite a com- 
plex dependence on the ratio w/T, yet the prefactor of the 
term vanishes at w = inT. We consider in detail 2D 
and 3D PL's, a 2D system at a nematic QCP, and also a 
marginal FL. We discuss under what conditions the local 
approximation is valid in all these cases. We also discuss 
in this Section how one should properly construct the 
self-energy along real frequency axis to make sure that a 
replacement of w by inT agrees with the analytical con- 



tinuation of the self-energy into the upper half-plane. We 
present our conclusions in Sec El 

In the subsequent papery^^ we discuss the constraints 
imposed by the first-Matsubara rule on the VI /T scaling 
the optical conductivity (j{Q.,T) of a FL, and the conse- 
quences of these constraints for the experiment. 

Throughout the paper, we denote the retarded self- 
energy along the real frequency axis as 'S^{uj,T) and 
the self-energy along the Matsubara axis as T,i^{ujjn,T), 
where k is the electron (quasi)momentum. We set the 
overall sign of the retarded self-energy via 



(1.1) 



where £k is the electron dispersion, and define the Mat- 
subara self-energy in such a way that it is real on the 
Fermi surface, i.e., 
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(1.2) 



II. SINGLE-PARTICLE SELF-ENERGY: 
CANONICAL FERMI LIQUID 

In this Section, we briefly review the derivation 
of the scaling forms for the self-energy in a conven- 
tional FL to order and T^: ImE^(a;,r) (x lo^ + 
TT^T^, Rei;^(a;,r) = Aw with no 0{u}T) terms in ei- 
ther of these quantities. We first show how these forms 
are obtained in the perturbation theory, then use the 
Eliashberg's argument^ to generalize the derivation to 
an arbitrary order in the interaction, and finally relate 
these forms to the first-Matsubara rule - a special prop- 
erty of the self-energy at the first fermionic Matsubara 
frequency ujjn=o,-i = ±7rT (Sec. Ill Cp . 

A. Perturbation theory 

We consider a system of fermions on a lattice with 
single-particle dispersion £k- We assume that the Fermi 
surface does not have nested parts and is away from the 
van Hove singularities but otherwise arbitrary. Near the 
FS, Ek can be approximated as ek = Vk^ -(k — k^), where 
ki? is a vector pointing in the direction of k and residing 
on the FS, and Vk — Vk£k- We will see that uj^ and 
terms in lm'E,^{uj,T) come from low-energies where the 
linear approximation is valid. Having this in mind, we 
follow a conventional reasoning of a FL theory, set the 
upper cutoff of the theory with the linearized dispersion 
at some energy A (generally comparable to the band- 
width, W) , and absorb all renormalizations from energies 
between A and W into non-singular renormalizations of 
the effective mass and quasiparticle residue Z. The bare 
Green's function of low energy fermions is then given by 

(uj) ^ Z^J {u - + i6) , (2.1) 



FIG. 1. (color on-line). Diagrams for the self-energy. K = 
(c.™,k), 7^' = (L^„,,k'), Q = (r!„,q). 



where (5 > is infinitesimally small and Z^^^, in general, 
varies along the FS. We further assume that fermion- 
fermion interaction, Uq, is static and non-singular for all 
q connecting points on the FS, including q — 0. This is 
the case for, e.g., a screened Coulomb interaction. 

The lowest-order diagrams which contribute to the 
imaginary part of the fermionic self-energy are shown 
in Fig. [TJ The imaginary part of the fermionic self- 
energy arises from the convolutions of two Green's func- 
tions marked by slanted dashes in Fig. [TJ In diagram 
a, such a convolution is just a particle-hole bubble, the 
imaginary part of which scales linearly with the bosonic 
frequency Q. In diagram b, this convolution involves 
the momentum-dependent interaction, but the result still 
scales linearly with 51. 

To see this in more detail, we write down a Matsub- 
ara form of the self-energy from diagram a) and obtain 
'E,^'k{uj,T) by analytic continuation. With our definition 
for the self-energy (|1.2p . we have 

S£ (c^„„ T) ^-iTY, [ C/(q)'Gk+q iuj,n + n,,) Hq (fin) , 

(2.2) 

where 

Hq(r!„)=2r^ f Gk' (w„,0 Gk'+q (c^,„' + f^«) (2.3) 
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with fl„ = 2TTnT and /i = / d^l/(27r)^. Performing 



analytic continuation in both (|2.2[) and (|2.3p . we obtain 
the retarded self-energy along the real frequency axis 



n + uj 



coth — G^Vq + r!) Im (Q) + tanh Im G^^q + K i^) 



2T 



(f]) ^2 ^[tanh ^ Im GJ; G^,+^ {io' + tanh ^^G^, {J) Im G^,+q + f7) 



(2.4a) 



.(2.4b) 



Extracting the imaginary parts of Eqs. ()2.4ap and (j2.4bp 
and using the relations tanh(a;/2) = 1 — 2np(x) and 
coth(x/2) = 2nB{x) + 1, where npix) and nsix) are the 



Fermi and Bose functions, correspondingly, we obtain 



ImS^^''(cj,T) 



xImG^+q(w + f7)Imn^(n) (2.5a) 
Imn^(f^) ^2 J J ^ [npioj' + npioj')] 

xImG^,(w')ImG^,+q(w' + 17). (2.5b) 
Equation (|2.5bp can be-rewritten as 



duj' [?if(w' + ri) — nF(w')] 



«k' (2^)2 



k' +q 



(2.6) 



where dAi./ is the element of the D — 1-dimensional 

FS. The integral over £k' gives 5{ijj' + — Ck'+q) |ek/=w' • 
The role of this (5-function is to impose a constraint on 
the angle between k' and q. Since this angle is not, 
in general, small, it suffices to resolve this constraint 
at oj' = ri = because, as subsequent integration will 
show, w' ~ ri '-^ max{aj,T}. The J-function thus re- 
duces to 5{e-\fj +q) |e^,=o , which means that both the ini- 
tial and final states are on the FS. [Notice that this ap- 
proximation corresponds to expanding the 5-functions in 
rasLyi{uj,T} / Ep rather than in w/ek-] The integral over 
w' now gives j'^^dbo' [npiyj') — npi^j^' + ^)\ = 5^, and 
Imllq (i7) reduces to f2 multiplied by a function of q, av- 
eraged over the FS: 



Imn^(f7) = - 



{2n) 



dA 



k' 



■^ki,-Z'k^+q<^(ek^+q) 



(2.7) 



For small q, the prefactor of fl behaves as 1/q. Substi- 
tuting Imllq (fJ) (X into (|2.5ap . and applying the same 
procedure as above to integrate over the momentum, we 
obtain, for k = k^? 



ImE«;"(L^,T) = 2G„ 



dnn [nB(ri) + nplui + fl)] 
(2.1 



with 



2(27r)^-ijqJ «k'. 



dA 



k' 



^kF+q^k'p^k^+q 



X(5(£kf+q)5(ek^+q)t/q- 



(2.9) 



[A factor of 2 in (|2.8p is introduced for future conve- 
nience.] The frequency integral in Eq. (12. 8p is readily 
evaluated 



dnn [nB{n) + np{uj + n)] = i (w^ + n^T^) 



hence 



ImE^'"(w,r) 



G^c 



2 



(2.10) 



(2.11) 



We can now specify what actually makes the analysis 
above applicable only to conventional FLs rather than to 
all FLs: it is an assumption that the integral in Eq. (|2.9p 
is convergent in the infrared. Power counting shows that 
the integrand behaves as l/g^ for g — > 0; the integral over 
d^~^q then converges for 13 > 2 and diverges for D < 2. 
Infrared divergence for D < 2 will modify the w and T 
dependencies of 'S^'kp{uj,T) compared to the canonical 
form valid for Z3 > 2. 

Diagram b is analyzed in a similar way with the only 
difference that the quantity [/qlmllq (il) in Eq. (|2.5ap is 
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replaced by 

ImP^k(f^)= / / ^K(c^' + r!)-nf.(e)]C/k_k' 

xlmG^,{e)lmG^,+^{e + n). (2.12) 

Still, Ini7^^k(^) scales as ft for — 0. Evaluating the 
integrals in the same way as above, we find 

lml]^;'(w, T) = Cb (w' + TT^T^) , (2.13) 

where 

X(5(£kf+q)5(eki.+q)[/qC/ki._k^. (2.14) 

As before, the integral in Eq. (|2.14p is convergent for D > 
2. Comparing Iml][^^°(a;, T) and ImS^.^ (aj,T), we see 
that they both have the same scaling form + 7r^r^ and 
differ only in prefactors which, in general case, depend on 
]i.F, i.e., on position along the FS. 

The real part of the self-energy can be obtained either 
directly, e.g., from Eq. (|2.4ap for diagram a, or via a 
Kramers-Kronig (KK) transformation of lmT,^{uj,T) 

2lu r°° JmE? (w',T) 
ReJ:ti^,T)^—V dLo' (2.15) 

where V stands for the principal part. The integral 
is ultraviolet divergent if Eqs. p. lip or p.l3p is used 
for ImEkp(a;,r), which implies that, to get the correct 
form of ReSkjj(u;,T) from the KK transformation, one 
has to use the full form of ImEky(a;,T) rather than its 
low-energy approximation. Nevertheless, one can easily 
make sure that, to quadratic order, ReSkj^ (i^, T) = \\s_pUj 
(where Akj? varies, in general, along the FS) with no loT 
term. 

A comment is in order here. By applying (|2.15p to 
(|2.1ip or (12.131) . we can only show that there is no "uni- 
versal", cutoff-independent luT term in ReYi^{uj,T)^ and 
hence no term in Rel]^(i7rr, T). There is stiU a pos- 
sibility that a term in ReE^(i7rr,T) may come from 
internal frequencies in (j2.15p comparable to the upper 
cutoff of the low-energy theory. We show later, in Sec. 
Ill C[ that this is not the case, and that only a term 
emerges from high energies. 

B. Arbitrary order in the interaction 

We now follow the argument by Eliashberg^ who 
showed that the o;^ -I- tt^T^ form of the self-energy at 
finite T holds to all orders in the interaction (a similar 
reasoning was also employed by Luttinger^- to show that 
ImE^kF(cj,T = 0) cxw^). 

The argument is as follows. In the second-order dia- 
grams, the -|- TT^T^ form comes from the region where 
all three intermediate fermions are located within the 



window of width of order w or T near the FS. Accord- 
ingly, the interactions C/q can be approximated by their 
values evaluated for the case when when the initial and 
final states are on the FS, i.e., q = Ij- — 1'^.. 

In a self-energy diagram of any order, one can select 
a cross-section with three low-energy fermions, and sum 
over all other fermions without assuming that they are 
near the FS. The diagrams of this kind can be cast in 
the form of Fig. [1]:). The three selected fermions are 
near the FS and the shaded squares are the full vertex 
functions. Because integration over the fermionic lines 
already gives a function quadratic in uj or T, one can 
set T = in the remainder of the diagram and project 
all four external momenta onto the FS. As long as the 
full vertex functions do not diverge, they do not affect 
integration over dispersions and frequencies of interme- 
diate fermions. Self-energy corrections to fermionic lines 
are also irrelevant because the dressed Green's function 
still has the form of Eq. (|2.ip at the lowest energies - 
adding one-loop self-energy to Eq. (|2.ip simply replaces 
iS by iC{iiP -|-7r^T^), which has an extra power of energy 
compared to uj and hence does not affect the and 
terms in the full self-energy. As a result, the + tt'^T^ 
form survives to an arbitrary order in the interaction - 
self-energy and vertex renormalizations only affect the 
overall factor in ImE^^(w,T). We then have for a con- 
ventional FL and to order aj^,r^ 

lmY.^^{uj,T) {uj^ + 7r^T^) . (2.16) 

The prefactor C depends on model parameters, including 
the cutoff A, and is thus non-universal. Substituting this 
form into KK formula, Eq. (|2.15p . and using the same 
arguments as in previous section, we find 

ReI]^^(a;,T) = Aw(l + X T) (2.17) 

(we spelled out the x T combination to emphasize that 
there KeT,^{u!,T) does not contain an uT term.) The 
prefactor A is again non-universal. 

C. Self-energy along the imaginary axis: 
the first-Matsubara-frequency rule 

We now show that the scaling form of ImE^^ (a;, T) in 
and the absence of the uT term in ReS]^^(w,r) 
are related to a particular behavior of the self-energy at 
the first fermionic Matsubara frequency LJm^o — ttT (the 
same behavior holds at ujm=-i = — ttT). 

1. Analytic continuation 

Let us first analytically continue ImE^(w,T) and 
Rei;^(w,r) in Eqs. ((2H1) and (pTfl) into the upper-half 
plane of the complex variable uj ^ z = z' + iz" . Because 
ImE^(aj,T) and ReI]^(a;,T) are analytic, their analytic 
continuation reduces to just a replacement of w by z. The 
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FIG. 2. (color on-line). Analytic structure of the function 
ImE {z, T) in the complex z plane. Analytic continuation 
from the real axis is possible to any point within the shaded 
region, including the points ±mT, but not beyond this region. 



complex function IniE^(z, T) cx + tt^T^ then vanishes 
dX z = inT and KeT,^(z,T) becomes inXT, so that the 
full self-energy reduces to iirTX and does not contain a 
term. This is only true, however, if z is replaced by 
first Matsubara frequency. For any other uj^n ttT, the 
function I]^{iu!m,T) contains a term. As we said in 
the Introduction, we will refer to this property as to the 
first-Matsubara- frequency rule (first-Matsubara rule). 

It is worth stressing that the rule formulated above 
applies only to the first Matsubara frequency. This may 
not seem to be the case if we replace w by in (2m + l)T 
with arbitrary m in ImS^^ (cj, T) given by Eq. (j2.5al) . 
before integrating over the bosonic frequency il. Doing 
so, and using the identity npi^ + iuJm) = —nsi^), we 
seemingly find that Im'S^{iu!m,T) vanishes not only at 
uj — ittT but also at any Matsubara frequency i'KT{2m + 
!)• 

This result is, however, false because the complex 
function, obtained by analytic continuation of, e.g., 
Eq. (I2.5ap . into the complex plane before the integral over 
n is performed, contains a sequence of branch cuts that 
run parallel to the real axis and intersect the imaginary 

I 



axis at the Matsubara frequencies (see Fig. [2]). As a re- 
sult, the imaginary part of the function 



F(z) = / dnn[nB{^) +nF{n + z)] 



(2.18) 



changes discontinuously a.t z = z' + i7r(2m -|- 1)T. For 
example, a discontinuity of Im_F'(z) at z = z' + inT is 



[Imi^(z' + iT(7r + (5/2)) 

= -6 dn 



sinh 



2 n-i 



2T 



lmF{z' + iT{7T~S/2))]\s^o 
2kTz'. (2.19) 



This implies that the substitution w — > z = z' -f iz" 
into the integral form of ImE^kF((u;,T), Eq. 12.81) . gives 
the same result for ImE^^(z,T) as the actual analytical 
continuation only in the region bounded by two branch 
cuts at z = z' + ikT and z = z' — inT, but not outside 
this region. In other words, the substitution w = iuJm 
into (|2.8p gives the correct result for only for the first, 
but not for all Matsubara frequencies. 



2. Direct proof of the first-Matsubara frequency rule 

The first-Matsubara rule can be also proven directly, 
by computing the self-energy for a conventional FL 
in Matsubara frequencies. For the electron-electron 
interaction, this was done in Refs. iandi; however, the 
proofs presented in these two papers are valid under two 
additional assumptions, namely, of small-angle scatter- 
ing and of a quadratic dispersion, ek = (fe^ — kp)/2m*, 
where m* is the renormalized effective mass. In fact, 
neither of these two assumptions are necessary. In what 
follows we first consider the case of arbitrary-angle 
scattering but still keep an assumption of a quadratic 
dispersion, and then generalize the argument for an 
arbitrary dispersion. 



Quadratic dispersion. To be specific, we consider 
the 3D case; other dimensions can be considered in a 
similar way. The clamshell self-energy diagram (diagram 
c in Fig. [T]) reads 



o x ^ x ^ f d^k' f d^p 
^^.A^,n,T) = -^Z T^T^u:„, j J^J 



(2.20) 



I 

where rk,p;k',p' is the renormalized vertex (a filled dia- to be near the FS; this assumption has already been used 

mond in Fig. [TJ:). Since we have already assumed that in Eq. (I2.20[) . 
the dispersion is isotropic, the Z factor is assumed to be 

isotropic as well. The momentum transfers can be arbi- To evaluate the momentum integrals, the dispersion 

trary, but all three intermediate momenta are assumed ep+k-k' needs to be expanded in Sp, Ek, and ek'- For a 
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quadratic dispersion, we obtain after some algebra 

9k, k 



Ep+k-k' 



2fcp . ^?k,k' 

— — sm I sm ■ 



cos Wp.k-k' 



+ep ^1 + 2 sin cos 6'p,k-k' 

, / , \ f n ■ 2 ^k.k' , . ^k,k' 

+ (Ck + £k' ) I 2 sm — ^ h sm — ^ cos Wp.k-k' | > 

(2.21) 

where 9i^m is the angle between momenta 1 and m. A 
similar analysis can be carried out for any isotropic but 
not necessarily quadratic dispersion. For small-angle 
scattering (?k,k' ^ 1, Eq. (|2.21l) reduces to a familiar 
form Cp+q = (p + vpqcosOp q with q = 2kp sin(6'k,k'/2)- 
For the momentum p + k — k' to be on the FS, the first 



term in Eq. (|2.21l) must be small; for generic values of 
^k.k'i this condition amounts to a geometric constraint 



cos f^p k-k' = — sm ■ 



7k, k' 



(2.22) 



or 0p.k-k' — iCTJ" + ^k,k')/2- We expand the first term in 
(|2.2ip around this value as 6'p^k-k' = ±(7r + 0k,k')/2 ~ oi 
with a ^ 1, and set a = in the remaining two terms. 
This gives 



ep+k-k' = wJ-fc_Fasin6'k,k' 



+£pCos 6'k,k' + (ek + ek') sin 



■ 2 ''k.k' 



(2.23) 



where v*p — kp/m*. Substituting the last result into 
(P^ . we obtain 



da- 



Vp 

1 



l{Q.n + W„i) - Ek' iWm' - ffp 



i{ujrrL' + ^n) - vphpusmOi^^i^' - Epcos^ 6'k,k' - (ek + £k') sin 



2 ''k k' kF,pj?;k'y ,pj?+kF-k^i k^,p_F+kF-k^;k_F,Pi? • 



(2.24) 



Constraint (|2.22l) is assumed to be imposed on the mo- 
menta entering both vertices in the last equation. The 
integral over a gives 



where dots stand for 0{T) terms. Summation over r2„ is 
straightforward, and we obtain 



v*pkF sin ^k k 



■sgn(u;„' -t- fin), 



(2.25) 



while the integral over Sp gives — iTrsgnWm'- Summing 
the product of the two sign functions over ujm' , we ob- 
tain a "local", i.e., integrated over the momentum, po- 
larization bubble as a sum of two terms: — |f2„|/7r and 
a constant, proportional to the ultraviolet cutoff of the 
theory. The constant contributes only to the 0{T) term 
in Ekp(7rT, T), and we consider it separately later. The 
|f2„| term is the one relevant to our purposes as we need 
to verify that Skj? (■"'T', T) does not contain a contri- 
bution. The prefactor of the |f2„| term is given by 



2C = 



m*Z 
27r 



«ok k' cos ——1 



xF 



k^,PF+kF-k'„;kF,pj; 



k_F,pF;k^,pF-l-kF-k^ 
(2.26) 



The remaining integral over £k' gives sgn(a;m -(- f2„), and 
the self-energy becomes 



SkF(Wm,r) 



-2C7rT^sgn(w„ 



(2.27) 



SkF(wm,T) = C(^2j,2_^2^) ^ _ _ (2.28) 
The term in the Matsubara self-energy obviously van- 
ishes for Wm = ±7rT. 



Arbitrary dispersion. Equation (|2.28p is also valid 
for an arbitrary fermionic dispersion, with the only differ- 
ence that the prefactor C now depends on the position 
on the FS. To see how this works, we expand Cp+k-k' 
near a FS. Knowing that £k and £k' drop out anyway, we 
set them to zero and expand tp+k-k' around pi? as 



£p+qjr £p-pF+p_F+q_F — £p_F+q_F + £f 



V. 



Jl 



, (2.29) 



where <\p = ^p ~ k^, uj' = vi • p^?, p^? = Tpp/pp, and 
we suppressed * m vp for brevity. 



Substituting this expansion into (|2.20p and replacing 
integrals over 3D momenta by integrals over the FS and 
over the electron energy, we obtain 
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1 



"Tkj? ,pp ;k'p .pi? +q_F Tk'y ,p_F +q_F ;k_F ,Pf I 



(2.30) 



where Z 



3 — 



PF+kF-k^ 



The condition that all 



three internal fermions are located near the FS implies 
that, at fixed qf, the angle between pp and must 
be such that the first term in Ep^+qp is small. Sup- 
pose that, at fixed q^?, the constraint £p0.i_|^q^ = is 
satisfied for a set of symmetry- related points on the FS, 
p^*. The vector pp spans a narrow solid angle around 
each of p^'; therefore we can expand the dispersion as 

EpF+qF = (Pf - Pf') • V, where v = VpO,>+qp. Since 
is still fixed, it can be chosen as the polar axis of 
a local spherical system, in which a point on the FS 
is described by an equation pp — r{0, 



are parametrized as pp — 



Vectors p^' 



Pp — r{9^'^ — a, 0°'* — /3) « p^p — are — Pr^ where rg and 

are the partial derivatives of r with respect to 9 and 
respectively, evaluated at the point (6''^'*, 0°'*). Suppose 
that V makes angle 7 with the polar axis and, without 
a loss of generality, assume that x axis belongs to the 
plane formed by vectors v and q^- Then, cos^pO.; ^ — 

cos 0°'* cos 7 and cos 6* oa ~ sa (cos 6'"'* -I- a sin 6*"'*) cos 7 

Pf ^ ' 

to linear order in a. Substituting all of the results above 
into Eq. (|2.29l) . we obtain 



correspondingly. 



a 



(p°/ sin 6l°''' - re cos 6*"'') cos 7 -1- Sp 



Pf 



(2.31) 

which generalizes Eq. (|2.23p for the arbitrary dispersion 
case. The measure of integration over the area dAp^ 
reduces to 

« ^Ppl sin e^'^dadp. (2.32) 



The rest of the calculations proceeds in the same way 
as for the quadratic-dispersion case; namely, integrating 
first over a, then £p and, finally, over ek, we reproduce 
the same product of the three sign factors as before. The 
final expression for the self-energy reduces to that in Eq. 
(I2.27P with a different prefactor, which varies over the 
FS. 

Interestingly, we found that the seemingly obvious 
result that frequency summation in Eq. ()2.27p yields 
(|2.28p can be reproduced only with a considerable effort 
if one uses the Euler-Maclaurin formula to sum over r2„. 
Namely, one has to keep not only the "conventional" 
terms with the integral over n and derivatives of the 
summand at n = 0, but also the "remainder" term which 



which is often neglected when the Euler-Maclaurin for- 
mula is applied in practice. We discuss this issue in 
Appendix |^ 



3. The linear-m-T term in the Matsubara self-energy 

Finally, we consider in more detail the 0{T) contri- 
bution to T,kp{TTT,T). For definiteness, we focus on the 
3D case and restrict to quadratic dispersion. If we inte- 
grate in Eq. (I2.20p over Sp and £k' in infinite limits, as 
we did earlier in this Section, and retain a constant term 
[denoted as 11(0))] instead of the \n\ term in the local 
polarization bubble, we obtain 



(2.33) 



where A oc 11(0). Because only the n = term con- 
tributes to the sum, Skp (ttT, T) = AT. 

This result holds only if we integrate over £p and £k' in 
infinite limits. Since, however, we have set the cutoff of 
our low-energy theory to A, integrations over £p and £k' 
should, strictly speaking, be performed between —A and 
A. The magnitude of A then depends on the ratio A/Ep 
and reduces to the previous result only for A ^ Ep. In 
the opposite limit of A ^ Ep, which is more appropriate 
for systems in which Ep is of the same order as the band- 
width, A is much smaller, namely, A ^ {A/Ep) ln{Ep/A). 
We show this in Appendix |Bl We also checked if there is 
a contribution to 'Eup {t^T, T) at finite A but found 
no such term. The next term after the AT is of order 
T'^/A^. This one is irrelevant to our purposes, as later in 
the text we show that in a generic 3D FL there are uni- 
versal terms of order {T^ / Ep)\n Ep /T, which are para- 
metrically larger than a non-universal term. 



III. SINGLE-PARTICLE SELF-ENERGY: 
NON-CANONICAL FERMI LIQUIDS AND 
HIGHER-ORDER TERMS IN CANONICAL 
FERMI LIQUIDS 

We remind the reader that the analysis in the previous 
Section relied on the assumption that the momentum in- 
tegrals, incorporated into the prefactor C in Eq. (|2.16[) 
for ImE^^(a;,T), are free from singularities. These in- 
tegrals include quasiparticle renormalization factors, the 
effective interaction between the quasiparticles, and the 
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prefactor of the Q term in the imaginary part of the po- 
larization operator [see Eqs. (|2.9p and (|2.14l) ]. The quasi- 
particle renormahzation factors and the effective interac- 
tion are non-singular at small q, but the prefactor of the 
n term scales as 1/g and may give rise to infra-red diver- 
gencies. The momentum integral in the expression for C 
is over the D — 1 components of q lying in a plane tan- 
gential to the _D-dimensional FS. This integral converges 
for D > 2, i.e., in a conventional FL, but diverges for 
D < 2, i.e., in a non-conventional FL. 

The issue we discuss in this Section is whether the first- 
Matsubara rule holds in a non-conventional FL, and in 
a conventional FL beyond the order. We will show 
in this Section that the next after the 0{T) term in 
i;k^.(7rT,r) scales as for any D, i.e., T.i,^{-nT,T) = 
XT + dT^ . 

The term is subleading to the one m a con- 
ventional FL (D > 2), and thus the first-Matsubara rule 
holds to order in this case. However, the leading terms 
in a non-conventional FL (1 < _D < 2) are also of the 
order, and thus the first-Matsubara rule does not hold 
in this case. In the next Section, we show that the first- 
Matsubara rule holds to all orders in T for any D near 
QCP, when the local approximation becomes valid. 

We consider first the marginal case of Z? = 2, and then 
discuss the cases of 2 < D < 3, D = 3, and 1 < D < 2. 



D = 2 



In D = 2, the self-energy is non- analytic: 
ImE^^(a;,r) cx w^lnlw] at T = and T^lnT at w = 
^]^g fjj-gt subleading term in ReE^ (w,r) 
scales as w|w| at T = and as T^sgnw for w ^ T}'^ To 
logarithmic accuracy, the scaling form of ImE^ (w, T) is 
given by^^^-^ 



ImS^^ (c^, T) = Ca (w^ + tt'^T^) In 



A 



(3.1) 



where C2 is a constant. By the KK relation. 



ReS« (^,T) 



Ao; 



sgnw (a;2 + tt^T^) . (3.2) 



At this level, the first-Matsubara rule is obviously sat- 
isfied. Beyond logarithmic accuracy, however, the situa- 
tion is different, as we will now see. 

Let us first calculate the self-energy in Matsubara fre- 
quencies. Consider diagram a) in Fig. [TJ The corre- 
sponding formula for the self-energy is given by Eq. (j2.2l) . 
We explore an earlier observatio n^^'^^" — that the non- 
analytic contributions to the fermionic self-energy come 



from forward- and backscattering rather than from scat- 
tering by an arbitrary angle. The internal structures of 
diagrams with forward scattering and backscattering are 
the same, i.e., it is sufficient to analyze only one of these 
two contributions. We consider forward scattering, i.e., 
focus on small momentum transfers g, and also assume 
that the FS is isotropic (a circle). Consequently, the self- 
energy does not depend on the position on the FS but 
we will still keep the subscript ki? which indicates that 
the self-energy is evaluated on the FS, as opposed to the 
self-energy evaluated away from the FS also considered 
in this Section. 

At small q, the polarization bubble behaves as 



nq(r!. 



(3.3) 



The constant term in 11 gives rise to an 0{T) term in 
SkF(j7rT, T). We neglect it for now but will re-instate it 
in the final result for Ekp(7rr, T). Keeping the dynamic 
part in p. 31) and introducing polar coordinates for mo- 
mentum integration, we obtain for the forward-scattering 
contribution of diagram a) to the self-energy at arbitrary 
momentum k 



Sk(Wm, T) = -iA2T V [ ^ 



qdqdcj) 



I a. I 



,(3.4) 



where A2 = 47rw^(0)/m and u{0) = TOJ7q=o/27r is the 
dimensionless coupling constant for forward scattering. 
Integrating over 0, we obtain 



I]k(Wm,T) = J 



dqq 
2tt 



sgn(cj™ 4- ri„) 



.(3.5) 



First, we discuss the self-energy on the FS. Substituting 
£k = into Eq. p.Sp and integrating over q up to A/vp, 
we obtain 

Skf (l^m,?") = 2" X! \^n\sgn(uJrn + ^n) 



X In 



For uj,„ — ttT, the last result reduces to 



(3.6) 



Sk,(7rr,r) 



mln 



2n-l/2 {A^+n^Y^^ + {A^ + {n + l/2'f) 



2n 



1/2 (A2 



,1/2 



(A2 



1/2)^ 



1/2- 

rr/2 



(3.7) 
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where A = K/{2ttT) ^1. To evaluate the frequency 
sum, we notice that the second fraction under the log- 
arithm is close to unity in both regions of m that are 
relevant for the sum, namely, for n <^ K and for m ^ A, 
when n ^ n^. In either case, 



(A2 + n2)'/' + (A2 + (n + 1/2)2) 



In 



1/2 



(A2 + n2)'/' + (A2 + {n- 1/2)2)^/' 2(A2 + n^) ' 

(3.8) 

With this simplification, the sum over n can be eval- 
uated exactly. Performing summation, and adding the 
0{T) contribution from the static part of the polariza- 
tion bubble, we obtain 



2'Kvj, 



K 



7rln2 



(3.9) 



where A ~ {A2h.)/v'^p is a non-universal constant and 
K = 0.9160 is the Catalan's constant (X -f 7rln2/4 ^ 
1.460). We see that I](7rr, T) does contain a universal, 
i.e., cutoff-independent, term. We recall that there is 
no such term in Z) > 2, when the self-energy is analytic 
to order . The presence of such a term in = 2 
implies that the first-Matsubara rule breaks down once 
the self-energy becomes non-analytic. 

For completeness, we also reproduced Eq. (|3.9I) by eval- 
uating first ImS^^ (cj, T) and then evaluating Skj,, (ttT, T) 
using the general KK relation between the Matsubara 
self-energy and ImS^ {uj,T) 



°° ImSf (uj,T) 

doj I ' (3.10) 







Applying spectral representation to Eq. (13. 4p and inte- 
grating over the momentum, we obtain for w > 



ImS;^^(^,T) 



A. 



X In 



VA2 - f72 + ^A2 + ujy 



uj\uj + 2Q\ 



(3.11) 



if w < 2A, and ImS^^ (w, T) = otherwise. To logarith- 
mic accuracy, this expression reduces to Eq. (j3.ip . 

Substituting p. lip into (|3.10p and setting — nT, 
we find that the main logarithmic term in ImE^ i^^T) 
contributes only to the 0{T) term in T,i,p{t:T,T). The 
violation of the first-Matsubara rule comes from the sub- 
leading 0(^2) and 0{T^) terms. We obtained the first 
term in p.9p analytically and reproduced the second 
term by integrating over uj in (j3.10p numerically. 

A complete expression for (ttT. T) to second or- 
der in the interaction contains contributions from dia- 
grams a) and b) in Fig. [TJ Each of these diagrams con- 
tains contributions from the interaction with momentum 
transfers equal to zero and to 2kp with amplitudes U{0) 
and U{2kF), correspondingly. Collecting all these contri- 
butions, we obtain a complete result for (ttT, T) to 



second order in the interaction asi^i^ 
Sk^(7rT,r) ^ttTA 

r2 



7rln2 



[3u2(0)-|-2u2(2fcF)-2u(0)u(2fcF)] Lft: , 

21l/ p y 4 

(3.12) 

where u{q) — mUt^/ {211). The combination of the cou- 
pling constants in (I3.12p can be expressed via the spin 
and charge components of the forward (/) and backscat- 
tering [h) amplitudes, F^'' and F^, defined by 

r^'';/35 = H^'^o.p5,s + r{'V„^ • (3.13) 

where subscripts c and s stand for "charge" and "spin" , 
respectively. To first order in [/q, 

F^ = -F{ = u(0), F;; = 2u(0) - u{2kF), F^ = -u{2kF). 

(3.14) 

Using these relations, one can re-express Eq. p.l2p as 



r2 



Ek.(^T,T)-^rA 
X [2{(F^)2 + 3(Ft)2} + (F{) 



K 



7rln2 
3(F{)2 



(3.15) 



Equation (|3.15p can be extended to a FL with an 
arbitrary interaction. One can show, using the same 
arguments as in Refs. [13 and IH, that the self-energy 
still contains the same combination of forward- and 
backscattering amplitudes, except for in a general case 
, and F;^^ are expressed not via u(0) and u{2kF) 
but rather via fully renormalized four-fermion ver- 
tices F(k, k; k, k), F(k, -k; k, -k), and F(k, -k, -k, k), 
which may depend on both transferred and total mo- 
menta. Explicitly, we have 



F^ 



f!: = 



F° 



2ti 

27r 
Z'^m 



F(k,k;k,k); F{ 



Z^m* 



27r 



-F(k,k;k,k) 



27r 



[2F(k, k, -k; k, -k) - F(k, -k; -k, k)] , 
F(k,-k;-k,k). (3.16) 



A complete expression for the self-energy at the first Mat- 
subara frequency is 



Ek^(7rT,r) =^^A- 
x[2{(F^)2 



T2 



K 



8Ef 
3(Ft)2}-K(F/)2 



7rln2 
3(F{)2 



m 
mZ 



(3.17) 



There is one additional complication: the result 
in Eq. (|3.17p is actually based on the expansion 
of the polarization bubble in frequency: for free 
fermions, this amounts to replacing (j3.3p by nq(il„) — 
— (m/7r) (1 — \V,n\/vFq)- The static part of nq(i7„) pro- 
duces the T term in p.l7p . while the (smaller) dynamic 
part produces the term. At weak coupling, one can 
safely set the lower limit of integration over q to zero, be- 
cause the contribution from the region q < CI/vf ~ T/vp 
produces only higher than terms. In a generic FL, 
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an expansion of the polarization bubble is possible for 
\rQn\/vFQ ^ 1, where F is the largest of the scattering 
amplitudes in (|3.15p . When F ^ 1, which happens either 
when the interaction is strong or when the system is near 
a Pomeranchuk instability,^^ the condition vpq ^ |Fr2„| 
sets a new lower cutoff for integration over q. We con- 
sider the case of large F in Sec. IIVI where we show that 
the existence of this cutoff affects the prefactor for the 
term in Eq. (|3.17l) . which gets smaller as F increases. 
The consequences of the first-Matsubara rule for the 
de Haas- van Alphen (dHvA) oscillations in a 2D FL were 
analyzed in Refs. and 0, where it was shown that the 
amplitude of these oscillations contains neither a In T 
nor a term resulting from the self-energy of quasipar- 
ticles. This result seems to contradict Eq. p.l7p which 
shows that the self-energy evaluated at LOm — t^T does 
have a term. In fact, there is no contradiction be- 
cause Eq. p.l7p refers to the self-energy evaluated on the 
FS, i.e., for £k = 0, while the dHvA amplitude contains 
the self-energy evaluated at the "Matsubara mass-shell" , 
defined by a solution of the equation G^^{uJm) =0. It 
turns out that these two self-energies do have different 
T dependencies. The amplitude of dHvA oscillations in 
any thermodynamic quantity contains the following di- 
mensionless combinatioiii^ 



A, 



dHvA 



iT 

2t:lOc 



(3.18) 

where oJc is the cyclotron frequency. For simplicity, 
we omit 0(ek) and 0{uJm) terms in Sk, which only 
renormalize the effective mass entering the cyclotron fre- 
quency, and focus on terms of order InT and higher. 
We also focus on the weak-coupling regime, when the 

I 



Matsubara mass-shell can be determined perturbatively; 
to lowest order in the interaction, the mass-shell simply 
coincides with the pole of the Matsubara Green's function 
Ek — i'^rn- Substituting the self-energy p.l5p evaluated 
at £k = i^m into Eq. (j3.18l) and integrating over ek, we 
obtain 



A. 



T 



dHvA 



^ exp -27r 



(3.19) 



where 



E(w™) = — I dqq2_^' 



2vr Jo 



(3.20) 



For high enough temperatures, i..e, for T > ujc, one needs 
to keep only the tUm = ttT term in the sum of Eq. (j3.19l) . 
This is where the first-Matsubara rule becomes useful 
because the Matsubara sum in Eq. p.20|) vanishes for 
u!m — ttT, and ^dHvA reduces to the free-electron result 
(modulo renormalized effective mass) with no extra T 
dependent terms. 

A related point is the difference in the behavior of the 
self-energy at finite T and at T = 0. At T = 0, the per- 
turbation theory in 2D for the self-energy diverges near 
the mass-shell, and needs to be resummed to eliminate 
these divergences .-i^iSS The mass-shell singularity shows 
up already in the second-order self-energy at T — 0, 
which is obtained by replacing the Matsubara sum in 
Eq. (|3.5I) by an integral over r2„. To logarithmic accu- 
racy, this yields 



Sk(a;„i,r = 0) = - 



A2 



+ - (w„i + iek) In 



A 



zek 



In- 



A 



- i£k 



(3.21) 



r 



The mass-shell singularity in this equation is manifested 
as a divergence of the first logarithmic term at £k = i^m- 
However, if we keep T in Eq. p.Sp finite, integrate over q 
at finite ek, and then re-arrange the resulting Matsubara 
sum, we obtain to logarithmic accuracy and for Matsub- 
ara frequencies with m = 0(1) 



A.T 



-ttT 



I]k(a;„i,T) = -^-^ ^ |f]„ 



(3.22) 



27rr 



X In 



A2 



{2Qn - i^m - i£k) [iOm + JEk) 



The limit of T — ?> in this equation reproduces Eq. p.2ip 
with the same mass-shell singularity at Ek = iiOm- How- 
ever, at iOrn = the sum in Eq. (I3.22p contains no 
terms and thus the mass-shell singularity in Sk(7rr, T) is 



absent. This is the reason why the mass-shell singularity 
does not show up in the dHvA amplitude. 



B. Higher-order terms in canonical Fermi Liquids 

(2 < D< 3) 

In canonical FLs, the first-Matsubara rule holds to or- 
der T^. Let us now verify whether if it also holds to 
higher orders in T. To obtain Ek(7rT,r) beyond the 
order, we need to go beyond the approximation we used 
in Sec. |TT1 where we assumed that the interaction con- 
nected only the points right on the FS. 

We verified that, as in 2D, the terms relevant to our 
analysis come both from small momentum transfers and 
momentum transfers near 2kp. Consider for definiteness 
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a small momentum contribution to diagram a) in Fig. [T] 
The corresponding formula for the self-energy is given by 
Eq. (|2?2|) . To single out potential terms in SkplvrTjT) 
beyond the order, we subtract from the integrand 
in Eq. (|2.2p its expression for the case when the effec- 
tive interaction connects the points right on the FS. We 
parametrize the measure of the D dimensional integral 
over q as d^~^q±dq^^, where a D — 1 dimensional vector 
q_L lies in the plane tangential to the FS and g|| is along 
the normal to the FS, and replace the integral over 
by that over the fermionic dispersion in the final state 
eicF-hq ~ vpq\\ = e. As before, we neglect the static part 
of nq(r2„), which contributes only to the 0{T) term in 
S, and approximate the dynamic part of nq(J7„) by the 

l^nl/? — I I / •y/'?! + 9f\ form. Using these simplifica- 
tions, we express the part of the self-energy not captured 



in Sec, nil as 

(5I](7rr,T) = 

1 



dq±de 



i2n) 



D 



(3.23) 

where Ad = 2i/£,7r-^/r[(£) - 1)/2]U^^„, T[x] is the 
Gamma-function, and ud is the density of states per 
spin projection in D dimensions. Because {t^T, T) 
obtained in Sec. Hi] contains only linear-in-T term and 
thus satisfies the first- Matsubara rule, potential devia- 
tions from this rule are due to 5Y.{nT,T). Integrating 
over q±_ and e in p.23|) . we find that 5Yj{'kT,T) contains 
a contribution 



T 



F n 

rpD 



= Ad' i2.f-^Y. 



(n+l/2)^-2g^ 



n+ 1/2 



-(n-l/2)^-2Q 



D 



n - 1/2 



(3.24) 



where A = A/2ttT, and 

Qd{z) (3.25) 
_^ dxx^-'^dy + + z2 _ ^a;2 _^ ^2 

Jo Jo (27r)^ {y^ + l)Vx^+7Jy/x^+y^ + ^' 

The sum in p.24p contains a contribution from the up- 
per limit, which just adds an extra piece to the 0{T) 
term, but it also contains a A-independent contribution 
from n = 0(1) which yields ST.{i7rT,T) cx T° . We see 
therefore that the fuU ^i,^{ttT,T) = 0(T) + 5J:{nT,T) 
contains a term, i.e., the first-Matsubara rule breaks 
down at order in conventional FLs. Still, Skj^ (ttT, T) 
and Y^it,p{uJm,T) for \ujm\ 7^ ttT are qualitatively differ- 
ent: the next term after T in Ekjj(a;m 7^ ttT, T) is 
while in Ekj,(7rT,r) it is T^, which for D > 2 is much 
smaller than T^. We verified that in the limit D ^ 2 the 
result matches the second term in (13. 9|) . For arbitrary 
2 < _D < 3, the sum has to be evaluated numerically. 

The case I? = 3 is special because Q3 diverges loga- 
rithmically. In this case we have, after integrating over 
X in (I3.25|) and neglecting non-logarithmic terms, 

5E(7rT,r) (3.26) 
= A,^Y.\^\{n + l/2) f\y^lyl±^, 

where z — n/{n+ 1/2). By power-counting, JE(7rT, T) 
scales as but there is an additional logarithm, which 
can captured by expanding the summand in Eq. (j3.26l) in 



1/n. The prefactor of the term is l/|n|, and the sum 
T^/\n\ yields a T^lnA/T term in Skp(7rT,T). Expand- 
ing the integrand of Eq. p.26p in 1/n, integrating over 
y, and collecting the prefactors for the T^lnA/T term, 
we obtain 

5I](^r,T) = A3-^ln^. (3.27) 

The complete result in 3D again contains the contribu- 
tions from diagrams a) and b) in Fig. [T] and includes terms 
coming from both forward- and backscattering. 

Note that the signs of (5S(7rT,T) are different in 2D 
and 3D [ cf. Eqs. dSJ]) and (13221)], i.e., the prefactor of 
the term vanishes at some D in between 2 and 3. 



C. Non-canonical Fermi liquids: 1 < D < 2 

The analysis for 1 < Z? < 2 parallels that in the 
previous Section. The extra term in the self-energy at 
LJm = ttT, given by (|3.24p . is still of order T^, and its 
prefactor is expressed via forward- and backscattering 
amplitudes. The only difference between the D < 2 and 
D > 2 cases is that, for D < 2, the term is larger 
than the one, and first-Matsubara rule breaks down 
completely, i.e., the next term after T in Skj^(wm,r) is 
of order for all ujm including = zLttT. 
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IV. THE FIRST-MATSUBARA-FREQUENCY 
RULE NEAR QUANTUM CRITICALITY 

A. Local approximation 

So far wc found that in a generic FL, either conven- 
tional or unconventional, the terms of order in the 
self-energy do not distinguish between the first and other 
Matsubara frequencies, i.e., the prefactor of the term 
in Y,\^p(LOm) is non-zero for all m. 

We now show that a different situation emerges when 
the system is tuned to the vicinity of a Pomeranchuk 
transition, at which a FL becomes unstable towards con- 
densation of particle-hole excitations with zero momen- 
tum transfer. A Pomeranchuk instability can occur in ei- 
ther the spin or charge channel. A magnetic (spin) insta- 
bility is likely to trigger pre-emptive transitions;^^^ and, 
to keep the discussion focused on the first-Matsubara 
rule, we only consider here a Pomeranchuk instability 
in the charge channel. In the bulk of this section we fo- 
cus on long-wavelength {q — 0) Pomeranchuk instability, 
(a quantum phase transition with dynamical exponent 
z = 3). At the end of this section, we briefly discuss the 
first-Matsubara rule near an instability at finite g in a 
system on lattice (a quantum phase transition with dy- 
namical exponent z — 2). 

Near a Pomeranchuk instability, interactions generate 
a large length scale ^ (the correlation length) which di- 
verges at the transition. In Z? < 3, a divergence in ^ 
brings the upper boundary of FL behavior down from 
0{Ep) to Cl>fl oc At large enough ^, cjfl becomes 

smaller than A, and the low-energy theory with the upper 
cutoff A describes now both the FL and non-FL regimes. 
We first consider the case of T ^ cjfl and then discuss 
the first-Matsubara rule at energies above wfl- 

The observation, which is most relevant to our anal- 
ysis, concerns the low-energy cutoff in the integration 
over bosonic momentum q in the formula for the self- 
energy, once we cast it into the form of Eq. p.4p . As 
we mentioned in Sec. IIII Al the term in Yi\^p{ujm) 
with a prefactor that does not show any special features 
at m = 0,-1 is obtained by setting the lower momen- 
tum cutoff to zero. This approximation can be justified 
at £, = Oil) ~ kp^, at least at weak coupling, but not 
at large ^. To show this, we follow earlier work^^ and 
assume that, near a Pomeranchuk instability with some 
angular momentum the fermionic self-energy given Eq. 
(13. 4p can be viewed as resulting from an exchange of low- 
energy and overdamped collective excitations. The prop- 
agator of these excitations at small q is given by 



Xo 



92 + e-2 + 7nq(f^„)' 



(4.1) 



where 7 depends on original fermion-fermion interaction 
and fermionic dispersion and, in general, is different for 
different As before, we keep only the dynamic part in 

nq(r!„). 



The one-loop self-energy is given by 

/d^q 
^^;^Gkf.+q(w™ + nn)xc^{^n)■ 

(4.2) 

An order-of-magnitude estimate for Yi\^p(uJrmT) can be 
obtained by expanding Xq(r2„) as 

Xq(^^n) = Xq(0) - X07f nq(f]„). (4.3) 

Substituting this expansion into Eq. (j4.2p and comparing 
the result to (|3.4I) in Sec. IIII Al and to its extension for 
an arbitrary interaction in Eq. p.lSp , we see that Xo7^^ 
plays the same role as the combination of the F^ terms 
in Eq. p.isp . i.e., the overall prefactor of the dynamic 
part of Yi\^p{uim,T) scales as 

Let us now look more carefully at the limits of in- 
tegration over g, which need to be imposed to ensure 
self-consistency expansion (|4.3p for Xq. Because nq(ri„) 
scales as |r2„|/5 at the smallest r2„, the expansion in 
nq(ri„) in (|4.3p holds only for q > 7|r2„|^^, which 
sets the lower limit in the integral over q. The up- 
per limit is set by Now, we expand the disper- 
sion as £kj7-i-q — vpq\\ + q\/2m* and express q\\ as 
% = S'kF+q./'^F — q±/i'2vFm*). Consider momentarily 
a free fermion propagator in (|4.2I) . Typical Ekj^+q are 
then of order Um + ^n, i-e., of order T for Matsubara 
indices m,n ^ 1. Since we expect the term to come 
from the region when both typical q\\ and q± are also 
proportional to T, the q']_/2m* term is of order and 
can be neglected compared to ek^+q- Hence typical g|| ^ 



£kF+q/vF, and typical q ^ y qi + [{uJm + ^n)/vF] ■ For 
large ^ and n,m ^ 1, (cj^ -I- Qn)/vF ^ ^n/vp is then 
smaller by than the lower limit for g, which is 7|f2„|f^. 
For 7t;_F^^ ^ 1, one can approximate q by q_L. This is 
equivalent to factorizing the momentum integral in (j4.2p 
as / dek^.-1-qGkf +q(wm + f^n) / d'^'^ q±x<iA^n) ■ In this 
approximation, the dynamic part of the self-energy in 
(1172)) reduces to 

Skf (t^™,T) = - — Vsgn(w,„ -|-n„)xL(f^n), (4.4) 
where 

XL(aO= / ^-^x^A^n). (4.5) 



Because XL{^n) is an even function of r2„, the r.h.s. of 
(|4.4p vanishes at |aj,„| = ttT, i.e., the first-Matsubara rule 
holds. For all other frequencies, such that |aj,„| ^ ttT 
but still \ujm\ ~ T, Ek^(cj™,r) behaves as T^^^u 
non-conventional FLs, and as 

Sk.(a.™,T) = T^e-^ {l + O {{ref-') } (4.6) 

in conventional FLs. 

We refer to an approximation, in which the momen- 
tum integral is factorized, as the "local approximation" . 
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The name reflects the fact that the fermionic self-energy 
in this approximation is a convolution of the density of 
states (the Green's function integrated over fermionic dis- 
persion) and the local susceptibility, obtained by integrat- 
ing the non-local susceptibility over D—1 components of 

If we keep uJm -I- i^n in g and compute Eky(a;m,7^) 
without making any approximations, we find that the 
term in I]]^^{ujjn,T) is present for all uJm', how- 
ever, its prefactor has different dependences on ^ for 

= ttT and all other uj„i. For \uJm\ ^ ttT, the 
term in Ti]s,p{ujm,T) is present even in the local approxi- 
mation, and the prefactor of this term scales as . For 
\^m\ = T^T , the prefactor is zero in the local approxima- 
tion, and scales as ^^(-0-2) j£ compute Skj? (ttT, T) in 
(|4.2p using a free fermion propagator. Using free-fermion 
propagator at large ^ is, however, not justified because 
the mass renormalization term, Awmj in is also pro- 
portional to ^ [this term involves a static susceptibility, 
Xl(0)]. Including this term into the Green's function af- 
fects the estimate for a typical ek^+q, which now becomes 
of order (1 -f )\)\ijJm + ^n\lvF- Accordingly, the prefactor 
of the T° term scales as A2^2(r'-2) |^^| ^ 
one-loop order, A cx ^-^--D (q^ In^ in = 3) and, hence, 
the self-energy at the first Matsubara frequency scales as 
Sk^ (ttT, T) cx T^^2 Still, for ah £> > 1, this is paramet- 
rically smaller than the self-energy at larger Matsubara 
frequencies, which, we remind, scales as T2^6-n 
ventional FLs, and as T^^^^ in non-conventional FLs. 

The main outcome of this analysis is that, near a 
Pomeranchuk instability the first-Matsubara rule ap- 
proximately holds, even if far from the instability this 
rule is broken, as it happens in non-conventional FLs. 
The distinction between the prefactors of terms in 
Skj?(wm,2^) likely persists to higher-orders in the loop 
expansion, even if higher-order corrections are not small. 
To verify this, we analyzed two and three-loop contribu- 
tions to the self-energy near a charge Pomeranchuk tran- 
sition in _D = 2. We recall that in = 2, the self-energy 
at a generic w,„ scales as T^f**. In 2D, a two- loop self- 
energy is small compared to the one- loop one, Eq. (|4.2|) . 
only if one extends the theory to N fermionic flavors and 
takes the ^ 1 limit J^iii The three- loop self-energy 
is not small even in the large- limit (Refs. [3 and 
[iTh , and higher-order terms even bring in additional log- 
arithmic singularities J^iiii^ We computed two-loop and 
three-loop contributions to the self-energy along Matsub- 
ara axis, and found that in both contributions the pref- 
actor for the term still vanishes at w„i = ±7rT if 
the local approximation is imposed, and scales as be- 
yond this approximation. Higher-order corrections may, 
in principle, generate additional logarithms and eventu- 
ally change the scaling of Skj? (""i T) with ^ from to 
with (3 < 2. However, because the one-loop results 
for Yj\^p{'kT,T) and Y,\^p(u}m,T) with |w„j| ^ ttT differ 
substantially (by a factor of in _D = 2), it is likely 
that the difference between the prefactors of terms 
in Ekjj (ttT, T) and in Ek^ (|i^m | 7^ ttT, T) holds to infinite 



order in the loop-expansion. 

The difference between Yi\^p{'kT,T) and Ekj:,(|wm| 7^ 
ttT, T) becomes particularly pronounced right at the 
Pomeranchuk instability. Now ujpL =0, and the self- 
energy exhibits a non-FL behavior at any finite lo or T . 
The self-energy for generic uJm 7^ ±7rT can be divided 
into two parts: dynamic, S'^, and static, Yi^. The dy- 
namic part comes from processes with non-zero energy 
transfers, corresponding to f2„ 7^ in the Matsubara 
sum of Eq. (|4.2[) . The critical form of the dynamic part 
is obtained by replacing in Eq. ()4.6|) by T, which 
gives E^^(a;„i,r) ^ j^d/z^ rpj^g static part comes from 
scattering of static critical fluctuations, corresponding 
to a single term with r2„ = in Eq. (|4.2|) . At finite 
^, this contribution behaves as E^^(r) cx T^^^^. At 
^ — 00, the static contribution diverges for D < 
This divergence is usually regularized by introducing a 
temperature-dependent correlation length, C(T), which 
remains finite at T > even right at criticality. On 
general grounds, one can postulate that ^(T) cx T^*^^ 
with Pt > 0, and hence E^^(r) cx t^-0t{3-d)^ 
the one- loop level, /3t = 1/2 (modulo logarithms) for 
2D quantum-critical systems with dynamical exponents 
Z = 2 and Z = but higher-order corrections may 
change the exponent. We will treat /St as a phenomeno- 
logical parameter of the theory. Comparing the expo- 
nents of the dynamic and static parts of the self-energy, 
we see that, for any D < 3, the leading T dependence of 
the self-energy is given by the dynamic part if /3t < 1/3 
and by the static part if /3t > 1/3. 

For the first Matsubara frequency, the static part of 
the self-energy is the same as for all other but the 
dynamic part is different. To obtain (ttT, T) at criti- 
cality, we re-evaluate the self-energy diagram in (|4.2|) by 
replacing the frequency in the denominator of the Green's 
function by the self-energy at the same frequency. Now 
typical £kF+q ^ Ekf(a;„i -I- r^„, T) |^^+n„^T = E(T). 
Expanding the bosonic propagator to leading (second) 
order in fk^+q and performing power-counting, we ob- 
tain 



E^ f7rr,r)«r- 



E cx 



^if < 1/3 

2;3T(3-D)^if > 1/3 ' 

(4.7) 



where we replaced E by E'^ and E*' for /3t < 1/3 and 
Pt > 1/3, correspondingly. We see that E^^(7rT,r) 
remains smaller than EJ^^ (wm ^ ttT, T) : the ratio of 
the two behaves as Y.^^{tiT,T) /T.'^ cx T'^(d~^)/3 for 
I3t < 1/3, and Ek^ (ttT, r)/Ed cx T<^D+i)/3~PTi3^D) 
Pt > 1/3. The exponent is positive for any Z? > 1 in the 
first expression and for 1/3 < (3t < {D + l)/3(3 — D) in 
the second one. In both cases the ratio of the self-energy 
at the first Matsubara frequency to that at a generic fre- 
quency scales to zero as T goes to zero. This smallness is 
a manifestation of the first-Matsubara rule at criticality. 
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Scaling form of the self-energy in the local 
approximation 



D = 2 



We again use (|4.ip for Xqi^n)- In D = 2 we have 



A non-trivial aspect of the first-Matsubara rule near 
a Pomeranchuk instability shows up when we consider 
the self-energy along the real frequency axis. At order 
, both the real and imaginary parts of Y,^^{lu,T) are 
rather complicated functions of uj and T, and the exten- 
sion of I]^^(cli,T — 0) to finite T by no means implies 
that UJ is replaced by \/uj'^ + tt'^T'^. Still, within the local 
approximation, we obtain, analytically continuing (|4.4p 
to real frequencies 



ImS« (c.,r) = 



2TrvF 



(4.8) 



ReE« 



The r.h.s. of Eq. (|4.8|) is an analytic function of complex 
variable uj ^ z = z' + iz" within the stripe |Imz| < ttT, 
see Fig. [21 Within this stripe, one can then analyti- 
cally continue ImE^^(a;,T) into the complex plane by 
just replacing uj ^ z. Aa a result, ImE^^ (ittT, T) 
is still given by (j4.8|) . but with zttT instead of uj in 
the r.h.s of this equation. Because nsi^) + ni?(il + 
mT) = 0, ImE^^(w = mT,T) vanishes. The full 
(zttT, T) vanishes by the first-Matsubara rule, hence 
uj,T) must also vanish (up to a 0(r) term), if 
we replace uj by ittT. These two requirements then set 
non-trivial constraints on the scaling functions of uj/T 
in ImE«^(w,r) oc |a;p//B(|c^|/r) and ReS^^(c^,r) cx 
a;|wp~VflD(|w|/T): both fioix) and fnoix) must van- 
ish at a; = itt. 

In the next two sections, we obtain explicit forms 
of //udwl/r) and fB.D{\^\/T) for near-critical FLs 
III D — 2 and D — i and show they they satisfy the 
constraint. 



Imxf(f^) - — — f71n— , 



(4.9) 



where wpL ^ 1/(7^'^) is the upper boundary of the FL 
behavior. The imaginary part of the self-energy is given 
by 



Iml]^„(a;,r) = Bq I rffifiln 



(4.10) 

where Bq = xol(,'^ /{'^'^'^Vp)- The real part of the self- 
energy is obtained via the KK relation. We skip the 
details of calculations and show only the final results. 
It turns out that the real part of the self-energy (the 
one which does not contain logarithms) can be computed 
exactly, up to the term of order uj which we omit below. 
The real part of the self-energy is an odd function of the 
frequency at k = ki;-. For w > we find 



ReS^^(L.,T) 

^_Bo 
4 



(4.11) 



t:uj'^ + inT'^ (—+U2 




)] 


\12 







where 



00 h. 



(4.12) 



fc=i 



is a polylogarithmic function. This expression can be cast 
into the scaling form ReI]J;^(w,r) = uj\uj\fR2{\uj\/T). 
For UJ = iirT, -e""/"^ = 1 and Li2(l) = 7rV6. 
Substituting these relations into (14. lip we find that 
ReS^^ (wT, T) = 0, as expected. 

The imaginary part of the self-energy is given by Eq. 
(|4.10p in the form of a one-dimensional integral. The 
formula for ImS^^ (w, T) can be simplified if we extract 
from it the leading logarithmic term. Combining the re- 
mainder of ImS^^ (w, T) with ReE^^ (w, T), we obtain 



^kA'^,T)^t—{uj +n T )\n—^ + i— [UJ 



4 

+2iBi 



In 



^0 



a;Li2(-e-"") 



1 



(4.13) 



where ln(— w^) = Inw^ — iTrsgno;. In Eq. (|4.13p . we sin- 
gled out the leading, logarithmic term in lmS^^{uj,T) 
and the rest has the form u)'^ fi2{\uj\/T). The scaling 
function is rather non-trivial, yet we see from (I4.13P that 
5]^^(i7rT, T) vanishes, as it should. 

In Appendix [Cl we discuss several subtle issues re- 



lated to analytic continuation of the self-energy to com- 
plex ui plane in a situation when either ReS^^(w,r) or 
lmS^^{uj,T) cannot be evaluated explicitly and has to 
be kept in an integral form, as in (|4.10p . 
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2. Subleading terms in D — 3 

A very similar situation emerges in 3D systems if we 
go beyond the leading, oj^ +7r^T^ term in the self-energy 
and consider the subleading terms of order and ut^. At 
T — 0, the real part of the self-energy scales as w'^ In 
and the imaginary part scales as At finite T, both 
parts contain scaling functions of |w|/r. The situation is 
somewhat similar to that in D = 2 in a sense that the 

I 



behavior is marginal due to logarithms. 
Using Xq(ri„) from Eq. (|4.ip we obtain 



8vF 



n\n\ 



(4.14) 



Substituting this form into Eq. (|4.8|) . we obtain after 
some algebra an explicit expression for ImS^^ (oj, T) to 



rder 



Im5]«^(c.,T) 



2rTi2\ 



Do 



N 
3 



TT^T^) +4T^ Li3(-e 



C(3) 



(4.15) 



where Co = xo7CV(327r4), Do = xo7CV(247r24), and 
C(x) is the zeta function. Using that Li3(l) — C(3), one 
can immediately verify that ImE^^ (iTrT, T) = 0, as it 
should. This happens despite that the functional form 
ImE^^(a;,r) is rather complicated at order uj^,T^; e.g., 
the prcfactor of the w"^ term is not the same as the pref- 
actor of the T"^ term. 

The real part of the self-energy contains logarithms 
and has to be left in an integral form. The calculation of 
Rei;^^(w,T) using the KK formula, Eq. (|2.15p . requires 

I 



some care as the integral is formally infrared divergent, 
if we use Eq. (|4.15p for ImI]^^(aj,T). The recipe is to 
i) start with the general expression for ImE^^(cij,T) in 
Eq. (|4.8p : ii) substitute it into the KK formula and ob- 
tain Rel]^^(w,r) in the form of a double integral; iii) 
keep the full form of Imx|J(f2) (without expanding it) 
at intermediate stages of the calculation, and change the 
order of integrations when it is convenient, iv) use the 
fact that 1uixl{^) vanishes in the infra-red and also 
that RexfW = (2/7r) dfflmxf (rj)/rj. Evaluating 
ReE^^(aj,T) this way, we obtain 



1 



2TrvF Stt'^vf 



xdx 



4a; 



Imx£(a;)da; 
ImXi (a;)x(ia; 



1 , X + UJ 1 
in ■ 



{x^ — uP'Y' 



2uj X — uj x 3(a;^ — o;^)^ 



ydy 



1 



e^/^ -I- 1 [(a;2 + y2 _ ^2-)2 „ 42,2y2 (2,2 „ ^2)2 



(4.16) 



where all integrals are to be understood as principal val- 
ues. The last two integrals are ultra-violet convergent 
for Imx^(x) given by (14.141) . The second term in (j4.16p 
is singular, but only logarithmically, and accounts for 



the w'^ln|w| term in ReS^^. Substituting Xl('^) from 
Eq. (pji)) and c ombin ing ImE^^ (w, T) and ReSk^ (w, T), 
we obtain from (|4.16p 



^Do 



Rexf(O) 
2tivf 
dx 



Do 



+ iCo (w^ + TT^T^) + —CO {uj^ + TT^T^) in 



37r 



X^ + Ul^ 



xix" -h TT^T^) -f 6T3 (Li3(-e-/^) - Li3(l)) 



(4.17) 



Equation (|4.17p is a complete expression for the self- One can easily make sure that I]^^(w,T) in (I4.17P is 

energy in a 3D FL within the local approximation. an analytic function of w in the upper half-plane, hence it 
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can be straightforwardly continued from the real axis into 
the upper half-plane just by replacing w by a complex z. 
KX, z — ittT, the second and third term vanish, while 
the last two terms cancel each other, i.e., at the first 
Matsubara frequency the self-energy contains a linear in 
T term but no terms of higher power of T, in agreement 
with the first- Matsubara rule. 



C. Marginal FL 

As another illustration, we consider the self-energy in 
a marginal FL (MFL)»22i The term marginal FL refers to 
a situation when the imaginary part of the self-energy 
is comparable to w, hence by the Landau criterion, the 
system is at the boundary between FLs and non-FLs. 
By the KK relation, if ImE(a;,r = 0) cx |a;|, then 
ReS(w,r = 0) (X wln(A/|a;|). 

Because in a generic non-conventional FL 
/mE^(a;,r) cx uj^ fioiH/T), the MFL behavior 
formally emerges when D approaches one. This limit 
is, however, special, and below we follow earlier work^^ 
and assume that the MFL behavior is associated with 
some sort of quantum criticality rather than with 
D = 1. Specifically, the MFL behavior emerges if 
one assumes Imx^(Sl,r) to be a scaling function of 
n/T such that Imxf(il,r = 0) = const x sgn(ri) and 
Imxf (f^,T) (X n/T for Q < T.^s A simple model form 
of Imx|^(ri,T) satisfying these conditions is 



flM M 



Imxf(n,T) = XLotanh 



(4.18) 



This expression is valid for ft smaller than some cut- 
off energy E* . At larger ft, lmxf^{fl,T) must decrease. 
To simplify calculations, we impose a hard cutoff, i.e., 
set Imxf (f^,T) to be given by for \n\ < E* and 

Imxf (f^,T) = for \n\ > E* . 

The first-Matsubara rule states that the self-energy at 
the first Matsubara frequency must be Ek^ (t^T, T) = 
ttTxl{0,T)/ {2ttvp). In all examples considered so far, 
we assumed that D > 1 and hence dropped this term, as 
it was of different order than the term which was our 
primary interest. Now ImE^^(w ^ T,T) = 0{T), and 
we should keep all 0{T) terms. 

Substituting Eq. (j4.18p into Eq. we obtain 



XLO 



dfl tanh — [nsifl] + 7ii?(Sl + uj)] 



2'kvf J T 

(4.19) 

Because the integral converges at large f2, and we are 
interested in a;,r ^ i?*, we can safely extend integra- 
tion over SI to the whole real axis. At T = we have 
from (|4.19|) ImS^^(w,0) = xlouj/27tvf, and at cj = 0, 
ImS^ (0,r) = XLQTTT/27TVF. When uj and T are both 




FIG. 3. (color on-line). Red: exact scaling function fiM{x) 
in ImE-^(a;,r) for the marginal-FL model. Eg. (|4:20]). Blu e: 
a scaling function obtained by replacing lj — >■ \/uj'^ + tv^T^ in 
lmS^{uj,T = 0). Square-root approximation (I4.21|l is practi- 
cally indistinguishable from exact fiM{x) in the interval of x 
shown in the figure. 



finite, integration in (14.19^ yields 

e^^ -I- 1 

Function // m (x) is plotted in Fig. [3l 

Expanding (|4.20p in w/T and casting the result into 
the form of a square-root, we obtain, approximately 



ImE«^(.;,r) 



XLO 
2'KVF 



,7r(4 - tt) 



(4.21) 



This form is obviously different from Vtt^T^ -|- o;^ ob- 
tained by by replacing w by Vtt^T^ -|- in the T = 
result. Nevertheless, substituting x = in into ()4.20p 
we find that fiM{™) vanishes, as it should by the first- 
Matsubara rule. 

The analysis of ReI]^^(a;,T) — ujfaMi^/T) requires 
more effort as one has to take care of the upper cutoff of 
the theory. The calculation is similar to the one we did 
for Z? = 3 in the previous Section. We use ImE^^(a;,r) 
in the form of Eq. (|4.19p , but keep the limits of the in- 
tegration over SI as —E* and E* and set E* to infinity 
only at the end of calculation. Without that, we would 
not reproduce the first-Matsubara rule for ReSkp (i^, 7")- 
Substituting lmT,^^{uj,T) from Eq. (|4l9)) into the KK 
formula we obtain after some algebra 
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fmiix) = - 



dy 



y 



ifiAiiy) -y) + - 



dy 



In 



y 



(4.22) 



where E* = E*/T. Both integrals are convergent and 
are easily evaluated numerically. 

At a; = in, the first integral yields (2/7r) x 0.96351, 
while the second integral gives {2/tt) (InE* + 1 — InTr) 
(2/7r) (lni5* — 0.14473), up to terms exponentially small 
in E* , which we neglect. Combining the two last expres- 
sions, we obtain 



fRM{i7T)^- {InE* + 0.81 



(4.23) 



According to the first-Matsubara rule, the result in 
Eq. (|4.23l) should be exactly the same as Rexf(0,r) 
[then Rei;^(i7rT,r) = i{T /2vF)fRM{iir) becomes equal 
to i(r/2wF)RexL(0,r)]. The static local susceptibil- 
ity is obtained by applying the KK transformation to 
Imxf (w,T) in Eq. 



I 

the higher-order contributions to I](cj„i, T) can be viewed 
as coming from processes with small momentum trans- 
fers, mediated by small q collective excitations of critical 
z — 2 modes, and higher-order contributions to S(a;m, T) 
from such processes do not vanish at iV = cx) (Refs. 
[H, [11 [H and [13). Still, in D = 2, a local propagator of 
the collective mode made of two z — 2 excitations scales 
as Xhi^n) oc J dqdq'dn',^x(q ,^'n)x{q + q' + ^'n) « 
ln|ri„| and is weaker than XLi^n) oc l/-\/|IXiI. As a 
result, the prefactor of the T^/^ term in J^{ujm,T), al- 
though does not vanish exactly at uJm — ivrT, is never- 
theless reduced by a factor of | InTj/r^/^. Contributions 
to this prefactor from even higher orders form series in 
|lnT|'7Ti/2 and may potentially give rise to an addi- 
tional anomalous power T''/Ti/2. The first-Matsubara 
rule then remains meaningful as long as rj < 1/2. 



Rexf(0,T) 



tanhx 



-dx 



JO 

2 



\nE* 



(In £;* +0.81878) 



dx \nx 
cosh^ X 

(4.24) 



again, up to terms exponentially small in E* . Comparing 
Eqs. (|4.23|) and (|4.24|) . we see that they are equal, as it 
should be, according to the first-Matsubara rule. 



D. Finite-g instability 

The discussion above is valid for a Pomeranchuk in- 
stability at q = 0, when the dynamical exponent z is 
equal to 3. In lattice systems, an instability may also 
occur at finite q, in which case z = 2, up to fluctuation 
corrections from multi-loop diagrams i^ii^i^ Such an in- 
stability is often called either spin-density-wave (SDW) 
or charge-density- wave (CDW), depending on whether it 
occurs in the spin or charge channel. The z = 2 case in 
more involved because typical q along the FS now scale as 
q± oc |fi„|-'^/^, while g|[ still scale as S(aJm + ^n,T). The 
one-loop self-energy for z = 2 problem scales as w^^^"^^/^, 
hence for i7„ ^ Um, typical q± are of order |r2„|'^^^^^/^. 
Local approximation is valid if typical g|| ^ q±, and is 
only justified for D > 2. At D — 2, q± and qn are of the 

same, order. The local approximation in this 

case can be imposed by extending the system to a large 
number of fermionic flavors N, and the analysis up to two 
loops indeed shows that the local approximation, and the 
first-Matsubara rule associated with it, become exact at 
N — oo. For a z = 2 transition, the first-Matsubara rule 
implies that I](w„i,r) evaluated at a generic Wm 7^ ±7rT 
contains a T^/^ term ( or a T^^'^ term in the FL regime), 
but the prefactor of this term vanishes at = ±7rT. 
The vanishing is not exact, however, because some of 



V. SUMMARY 

In this paper we analyzed in detail the fermionic self- 
energy 'E{uj,T) in a FL at finite temperature T and fre- 
quency uj. Our main goal was to understand how gen- 
eral is a certain property of the self-energy, the first- 
Matsubara-frequency rule. This rule states that the self- 
energy S(a;m, T), evaluated at discrete Matsubara points 
i^m = 7rT(2m + 1), exhibits a special behavior at the 
first fermionic Matsubara frequency namely, S(7rr, T) 
does not contain terms higher than OT. As a partic- 
ular manifestation of this rule, the imaginary part of the 
self-energy on the FS in a conventional 3D FL behaves 
as lmT,{uj,T) cx + tt^T^, with exactly a tt^ factor 
in front of the term, and ReS(a;,T) contains an lu 
term but no wT term. We found that the rule is not 
an exact one, i.e., Yi{'kT,T) in a generic FL does con- 
tain higher than linear terms in T . Still, the first term 
after 0{T) in 1:{'kT,T) in any dimension 1 < D < 3 
is of order (T^lnT in 3D). In D > 2, this term is 
parametrically smaller than term which is present in 
E(a;™,T) for \ujm\ ttT. We found that the term 
comes from only forward- and backward scattering, and 
is expressed in terms of fully renormalized amplitudes 
for these processes. We further showed that the first- 
Matsubara-frequency rule becomes exact in the local ap- 
proximation, when the interaction can be approximated 
by its value for the initial and final fermionic states right 
on the Fermi surface. In this approximation, which is 
justified, e.g., near a Pomeranchuk instability even if the 
vertex corrections are non-negligible, the term and 
all higher order terms in T,{'kT,T) vanish, and only the 
0{T) term survives. The first-Matsubara-frequency rule 
then imposes two constraints on the scaling form of the 
self-energy: upon replacing w by inT, lmY^^{Ld,T) must 
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vanish and ReS^(a;,T) must reduce to an 0{T) form. 
We considered several examples of the first-Matsubara 
rule, and argued that these two constraints should be 
taken into consideration in extracting scaling forms of 
E^(a;,T) from experimental and numerical data. 
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Appendix A: Evaluation of the Matsubara 
self-energy using the Euler-Maclaurin summation 
formula 



In this Appendix, we show how to reproduce the first- 
Matsubara rule for the fcrmionic self-energy by using the 
Euler-Maclaurin (EM) formula for summation over the 
Matsubara frequencies. Unexpectedly, the calculations 
involving the EM formula turn out to be quite involved, 
and to reproduce the first-Matsubara rule one has to keep 
not only the "conventional" terms in the EM formula, 
with the integral over a bosonic Matsubara frequency 
fin and the sum over the derivatives of the summand 
at n = 0, but also the remainder term, which is often 
neglected when the EM formula is applied in practice. 

To be specific, we consider Eq. (|4.2I) for the self-energy 

'^kpi'-^m, T) and set = ttT, which gives 



iTTT(2n + 1) — uf^ii 



(Al) 



with q — (q||,qi^). We assume that the local approxima- 
tion is valid, i.e., that typical are small compared to 
typical g_L, and the dependence of the bosonic propaga- 
tor on gii can be neglected. Within this approximation. 



Eq. (jAll) simplifies to 

Sk,(^r,T) = -^XL(o) 

T 



(A2) 



— ^ XL {2nT72) [sgn(2n + 1) - sgn(2n - 1) 



n=l 



where 



XLi27rTn)=. j d''-\^XciA'l^nT) / {2n] 



D-l 



(A3) 



The first term in (|A3|) is proportional to T, and the sec- 
ond term vanishes identically because for any n > 1, 
sgn(2n +l) = sgn(2n - 1) = 1. Hence, Sk^(7rr,r) does 
not contain terms beyond ©(T), in accordance with the 
first-Matsubara rule. 

An unexpected complication arises when one attempts 
to reproduce the vanishing of the second term in (jA2p by 
applying the EM formula to the sum over n. Under the 
condition that x(x) and its derivatives vanish at a: — ^ oo, 
which we assume to hold in our case, the EM formula 
reads^ 



/.oo 

1 JQ 



/(O) 



N 



2p- 



0=1 



/(2p-l)(0) 



Rn 



(A4) 



where Bk are the Bernoulli coefficients, /*^"-' is the n-th 
derivative of /, and Rm is the Poisson remainder term 



Rn — 



f 

Jo 



Bn{{1 



Nl 



(A5) 



where B^^x) is Bernoulli polynomial, and {1— .t} denotes 
the fractional part of 1 — x. 

In applications of this formula, it is often assumed that 
the remainder term Rj^ tends to zero in the limit N — 
oo and is thus dropped. We show that in our case the 
remainder term cannot be neglected and one should use 
the full EM formula, Eqs. (R4l and \Mil instead of the 
truncated one. 

Indeed, in our case, 

fix) = ^Xl{2tiTx) [sgn(2a; + 1) - sgn(2.x - 1)] . 

2VF 

(A6) 

The first term in the r.h.s. of (jA4p is the integral 
/o f{x)dx. Integrating f{x) from (|A6P over x we ob- 
tain 

1/2 



rOO rj-i /'1/2 

/ f{x)dx = — / dxxhi'^irTx) 
Jo Jo 



(A7) 
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where dots stand for the terms of higher order in T. 
Combining <\A7\ with the boundary term — /(0)/2 = 
— {T/2vf)xl{^), we see that the Unear-in-T term can- 
cels but the quadratic term contributes 



Si^tVlWt— (A8) 



to i]kr(7rr,r). 

This term would violate the first-Matsubara rule 
and must be canceled by the terms with the derivatives 
/("^(O). Because the derivatives of sgn(2a:: + 1) and of 
sgn(2a; — 1) vanish at a; = 0, one has to differentiate only 
XL (^ttTx). a contribution to Ek^ (ttT, T) comes from 
the first derivative of f{x), i.e., from the p = 1 term in 
the sum over p in the r.h.s. of (jA4l) . Terms with p > 1 
contribute higher powers of T. Using that B2 — 1/6, we 
find the contribution from the infinite sum with the 
derivatives /'■"■' (0) as 



S2 = -TVl(0)7^ (A9) 



and applying it TV — 2 times, we rewrite Rn as 

Rn = ^XLiO) [ B2[l-x)^ [x{l - sign(2a; - 1))] 
2vF Jo dx^ 

(All) 

Using that d/dx [I - sgn(2a; - 1)] = -25{x - 1/2) and 
also that Bi(l/2) ^ an d ^2(1/2) = -1/12, we obtain 
after integrating in (jAlip by parts 

S3 = -Rn = — x'l(0)S2(1/2) = -TVl(O)-^. 

Vf IZvf 

(A12) 

Combining the three contributions, we see that S1 + E2 + 
E3 — 0, as it should. 

An alternative way to compute the sum over bosonic 
Matsubara frequencies using the EM formula would be 
to "smear" the discontinuity in f{x) by integrating over 
<Z|| in (|A1[) in finite limits —Q < q\\ < Q and take the 
limit Q — > 00 only at the last stage. In this scheme, the 
remainder term i?oo does not contribute, but terms with 
p ~ ttT/Q become relevant in the sum over p in (jA4|) . 
This calculation is, however, more involved than the one 
we presented above, and we did not find a clear proof that 
the contribution from p ^ ttT/Q exactly cancels Ei + S2. 



The sum Ei + E2 = r2^^(o)7r/(12wi^) is non-zero. 

The sum over p in (jA4[) can be safely extended to in- 
finity as only the p = 1 term contributes a in the 
self-energy. If we used the truncated EM formula with- 
out the remainder term, we would have then obtained an 
incorrect result that Ek^ (ttT, T) does contain a term. 
In fact, the counter-term canceling the parasitic Ei -I- E2 
contribution does come from the remainder term in 
(|A5|) . even if we take N = 00 limit. Indeed, let's focus 
on term in the self-energy and replace xl{'^'^Tx) by 
2ttTx'l{0)x. We then have 



r2 

fix) ^ TT—x [1 - sign(2a; - 1)] . (AlO) 
2vf 

One can easily make sure that f{x) and its derivatives 
are non-zero only in the interval < x < 1/2, where the 
fractional part of 1 — x in the argument of the Bernoulli 
polynomial in (lASp is equal to just 1 — x. Furthermore, 
the derivatives /^"■'(x) with n> 2 vanish at the bound- 
aries of the integral in (jA5|) . hence one can integrating 
by parts N ~2 times and the boundary terms. Using the 
property of Bernoulli polynomials B'^{x) = N Bn^i(x) 



Appendix B: The dependence of the self-energy on 
the upper cutoff of low-energy theory 

In this Appendix we show that the prefactor of the 
linear-in-T term in the fermionic self-energy at the first 
Matsubara frequency Ekj^ (ttT, T) = AT depends on the 
ratio of the Fermi energy Ef — VFkF/2 to the up- 
per cutoff of the low-energy theory denoted as A. The 
result shown in Eq. (|2.33p with the "mass renormal- 
ization factor" A oc 11(0) corresponds to the situation 
of A <C Ep, when integration over intermediate ener- 
gies in the expression for the self-energy, Eq. (|2.30l) . 
can be extended to infinity. In the opposite limit of 
A <^ Ep, X is much smaller. To see this, we note that 
typical momentum transfers q = |k — k'| are of order 
kp, hence typical internal energies in the self-energy dia- 
gram are of order Ep- The integration over e^' in (|2.30p 
in finite limits changes the factor of sgn(a;m + ^n) to 
(2/7r) arctan [A/(ajm -f f2„)], which becomes small when 
typical ilm ^ Ep is much larger than A. The polariza- 
tion operator also changes, but the 11(0) term remains 
the same because it comes from the smallest frequen- 
cies. To simplify the computations, we keep nq(r2„) in 
the same form as before, but replace vpq by Ep, i.e., we 
set Uqi^n) = -n(0)(l - \nn\/y/nl + Ej,). Substituting 
this expression along with the result of integration over 
Ek' into the self-energy, we obtain 
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|»n| 



arctan 



A 



= XT- 



arctan 2 A + ^ I 1 



n = l 



A A 

arctan — — arctan 



71+ 1/2 



71-1/2 



(Bl) 



where A cx jllol is the same as in (|2.33p . A — K/2ttT 
and Ef = Ef/2ttT. For A :$> Ep, the term (1 - 



n/y'n? + Ep) decreases rapidly for n > Ep, when the 

difference between two arctangent functions is stih small, 
of order Ep /A. Then the first term in the last line in (jBip 
is the dominant one, and using that arctan 2A « 7r/2 one 
recovers Ekp_(7rr,T) = AT with A = n(0). In the oppo- 
site limit oi Ep ^ A, the first term can be approximated 
by unity for all 7i up to n ~ Ep ^ A. Because the dif- 
ference of the two arctangents scales as I/tt,^ for n ^ A, 
and 



E 

n=l 



arctan — 



A 



1/2 



arctan ■ 



A 



1/2 



arctan 2A, 



(B2) 



the contribution to the sum from positive n almost 
cancels that from n = 0. A straightforward analy- 
sis shows that in this limit, A is small and scales as 
A - {A/EF)\nEp/A. 



Appendix C: Analytic continuation of the 
self-energy in the local approximation in 2D 

In this appendix we discuss some subtleties of analytic 
continuation of the self-energy into the complex uj plane 
in a situation when either ReE^(aj,T) or ImY,^ {uj,T) 
cannot be evaluated explicitly and have to be kept in 
an integral form, as in (|4.10p . One can use the fact 
that the integral converges in the ultraviolet and mod- 
ify the integrand by shifting the variable. By doing so 
one can obtain several different formulas forImI]^(w, T), 
which all nevertheless yield the same result along fre- 
quency axis. The danger of this trick is that, by shifting 
variables, one imposes the dependence on the external uj 
onto Imxl^ which, in < 2, is a non- analytic function 
of its argument. As a consequence, if one now performs 
analytic continuation just by replacing w — >■ z, one ob- 
tains a branch cut which stretches down tolmz — >■ 0, and 
the self-energy will not obey the first-Matsubara rule at 
z = inT. To To make sure that this rule is satisfied, one 
has to use the Cauchy formula for analytical continuation 
which, in this case, is not equivalent to just replacing lo 
by a complex z. As an illustration, we consider Eq. (|4.10p 
in the form it was presented in Ref. [29l 



IniE^(a;,r) 




eA2 



uj^ In- 



I.UITtt'^T'^ + 2 



UJ In 



x — CJ 



X + UJ 



a; In ■ 



(CI) 



Along real frequency axis, this formula yields exactly the 
same result as Eq. (|4.13l) . However, if we formally re- 
place UJ by inT + S, with infinitesimally small S > 0, 
before integrating over x in the last term in (jCip . we 
obtain (i?o7r^T^/4)(0.17 + in) which obviously does not 
satisfy the first-Matsubara rule. The reason is that the 
w dependence is under the logarithm in the last, inte- 
gral term of (|C1|) , and each of the two logarithms there 
has a branch cut. Let's set uj = nTe^'^ and vary be- 
tween zero ( the real axis) and 7r/2 (the first Matsub- 
ara frequency along the imaginary axis). To understand 
what is going on, it is enough to move only little off 
the real axis, i.e., to consider only small (j). The branch 
cuts in the first and second logarithms in the integral 



term of Eq. (jCip are sA, x — \ and x — cos2(/), corre- 
spondingly. Each of them gives rise to a discontinuity 
in the imaginary part (Inz — In |z| -I- zTrargz, and the ar- 
gument of z changes discontinuously at the branch cut). 
At (/) = 0, the discontinuities coming from the two loga- 
rithms cancel each other, but at finite there is a range 
of X in between 1 — 0^/4 and 1, where the arguments 
add up to almost 27r This additional contribution makes 
ImE^(T, TrTe*"^) to be different from the one obtained by 
analytical continuation of Eq. (jCip . For small 0, the dif- 
ference is (tt^TVS) X {2m){(j)'^ / 4) / [e^ + 1). We verified 
numerically that this expression is exactly the difference 
between the analytical continuation of (|C1I) and the brute 
force replacement ttT -j> TrTe*"^ in ((CT|) . 
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